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We propose a new procedure for constructing projected three-pion correlation functions which 
reduces undesirable artificial momentum dependences resulting from the commonly used procedure 
and facilitates comparison of three-pion correlation data with theoretical models. 
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I. INTRODUCTION 



Two-particle intensity interferometry, exploiting Bose- 
Einstein correlations between pairs of identical bosons 
(pions, kaons), has been extensively used to extract in- 
formation about the space-time structure of high-energy 
hadron-hadron and nucleus- nucleus collisions . Three- 
particle Bose-Einstein correlations were shown to yield 
additional information which can not be extracted from 
two-particle correlations 0,0,0,13: (i) After subtracting 
the two-particle correlation contributions and extrapolat- 
ing to zero relative momentum, the strength of the gen- 
uine three-particle correlation term relative to the two- 
particle correlator provides an unambiguous measure for 
the degree of chaoticity of the source (or, conversely, it 
limits the degree of phase-coherence in the particle emis- 
sion process) This idea has recently been applied 
by the STAR Collaboration ^] to three-pion correlations 
from 130 A GeV Au-I-Au collisions at RHIC, showing that 
at these collision energies the data are consistent with a 
completely incoherent ("chaotic") pion emission process, 
(ii) For a completely chaotic source, the full momen- 
tum dependence of the "reduced" three-pion correlator 
r3(pi,p2,P3) can be used to extract spatial asymmetries 
of the source around its point of maximum emissivity 
0, for a partially coherent source, the momentum 
dependence of r^ constrains the relative spatial sizes of 
the coherent and incoherent parts of the emission func- 
tion 0. This information is, however, difficult to ex- 
tract because the required accurate analysis of its full 
momentum dependence puts statistical demands on the 
measured three-pion correlator which can not be met by 
presently available data. 

Statistical limitations therefore so far force experi- 
mentalists to project the three-pion correlation function 
onto a single relative momentum variable. The preferred 
choice 0, la, 13, for this variable is a Lorentz scalar that 
is completely symmetric under interchange of the three 
pions: Qs ^ y/ QI2+Q23+QI1 where Qf^ ^-{p^-pj^ are 
the Lorentz-invariant relative momenta between pairs in 



the pion triplet. This projection has undesirable, but 
to some extent unavoidable consequences for the rela- 
tive momentum dependence of the three-pion correla- 
tor. Specifically, it is known 0| that for a completely 
chaotic source the leading relative momentum depen- 
dence of r3(if, qi2, 923) (where K = (pi+p2+P3)/3 is 
the average momentum of the triplet while qi2 =Pi—p2 
and <723 =P2— P3 are two linearly independent relative 
momenta) at small is of 4*'' order in the compo- 
nents of the relative momenta . The projection on 
instead introduces a dominant quadratic g-dependence 
which buries the coefficients of the 4*'^- and higher-order 
contributions that would allow to extract additional in- 
teresting source information, such as the momentum de- 
pendence of the point of maximum emissivity and the 
asymmetry of the source around that point |3| . 

In this short note we investigate in some detail the ef- 
fects of this projection procedure. We show that there 
are different possibilities to project onto a single vari- 
able Qa some of which exacerbate the above problem far 
beyond the unavoidable minimum. In particular, the ver- 
sion employed so far in the data analysis has the undesir- 
able feature of introducing a relative momentum depen- 
dence in the projected three-pion correlator even for a 
static Gaussian source for which the underlying unpro- 
jected correlator is completely momentum-independent. 
We therefore suggest a different projection method which 
avoids this undesirable feature and furthermore allows for 
a "staged projection" onto any subset or different com- 
bination of the three particle momenta. As a side effect, 
the new method also facilitates the comparison of the 
projected three-pion correlator with theoretical models. 



II. TWO- AND THREE-PARTICLE 
CORRELATORS 

The two- and three-particle correlation functions are 
defined as the ratios between the two- and three-particle 
coincidence cross sections and the products of indepen- 
dent single-particle cross sections: 

^2(Pl,P2) 
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C2(12) = C2(P1,P2) = 
C3(123) = C3(pi,i>2,i>3) 



iVi(pi)7Vi(p2) ' 

^^3(Pl,P2,P3) 
^l(Pl)A^l(P2)A^l(P3) 



(1) 



2 



where 



N2{Pl,P2) 
^3(Pl,P2,P3) 



(Z = 1,2,3), 



(Ppi(Pp2 ' 
(Ppi(Pp2<Pp^ 



(2) 



Equations |^ are frame-independent even though they 
don't look that way since the energy factors Ei mak- 



ing the cross sections EidN / (Ppi etc. Lorentz invariant 
cancel between numerators and denominators. In the ab- 
sence of correlations C2 = 1. One defines the "true 
n-particle correlators" i?„ (cumulants) by subtracting 
uncorrelated parts and correlations involving fewer than 
n particles: i?2(ij) = CzfeO-l and i?3(123) = C3(123) - 
i?2(12)-i?2(23)-i?2(31)-l. The "normalized true three- 
particle correlator" r3(123) =r3(pi,p2,P3) 111 is ob- 
tained by dividing i?3(123) by the square root of the 
product of 2-particle correlators: 



^^(^23) ^ (g3(123)-l) - (C2(12)-l) - (g2(23)-l) - (C2(31)-l) 
y^(C2(12)-l)(C2(23)-l)(C2(31)-l) 



(3) 



r 



For a fully chaotic source, approaches the value 2 at 
zero relative momenta, qi2 = 923 — 0, for any value of the 
triplet momentum K = (pi+P2+P3)/3. The leading de- 
viations from this limit in both numerator and denomi- 
nator are of second order in the components of , but 
these leading corrections cancel in the ratio, leaving a 
leading 4'^ order term in rs 0| . It is obvious that a pro- 
jection of the numerator and denominator onto a smaller 
number of momentum components will in general destroy 
the cancellation of the leading 2"*^ order terms, thus gen- 
erating an artificial leading quadratic q^j -dependence of 
the projected r^. This leading quadratic g-dependence is 
not interesting because its coefficients depend on prop- 
erties of the source which can already be extracted from 
the two-particle correlation function. Unfortunately, this 
extraction is not completely model-independent , so it 
is not clear whether, even in principle, it might be possi- 
ble to subtract the unwanted second order terms, using 
the measured two-particle correlators, in order to uncover 
the subleading 4*^ order terms in the projected r^. We 



will not address this issue here. Instead we will focus 
on a different question: Is it possible to minimize these 
projection artefacts by using an optimized projection al- 
gorithm? 

III. SHORTCOMINGS OF THE "STANDARD" 
PROJECTION METHOD 



In |g| the STAR collaboration constructed a projected 
version f^^Q^) of the normalized true three-pion correla- 
tor using the following procedure: First, they constructed 
projected two- and three-particle correlation functions, 
C2{Qij) and C^lQs), by taking pairs and triplets of pions 
from the same collision event, summed over all events, 
binning them in the Lorentz-invariant variables Qij and 
Q3, respectively, and dividing the contents of each bin 
by that of an equivalent bin at the same Q-value which 
was filled with uncorrelated pairs or triplets generated 
from mixed events ||6|. This procedure can be formally 
represented by the following definitions: 



C2iQ.j) = 

C'3(Q3) = 



/ d^p.d^pj N2{p^,pj) 5 (Qy - y^-{pi-Pj)A 

S \ ' = 1,2,3) 

/ #p,d3p, m{p,)N,{p,) S [Q,, - v/-fe-p,)2) 

/ d^pid^p2d^p3 N3{pi,p2,P3) S (Qs - \J-{P\-P2f - {P2-PzY - (P3-Pl)^) 

/ d^pid^p2d^P3 Ni{pi)Ni{p2)Ni{p3) S (Q3 - \J -{pi-p2f ~ {p2~-pzf - (P3-Pi)^) 



(4) 
(5) 



From these projected correlation functions one then constructs the ratio ||2J) by fixing Q3 and summing over all 
two-particle relative momenta whose squares add up to Q\. Formally 



^3(03) 



^ dQi2dQ23dQ3i S(^Q3 - \J QI2+QI3+QI1 



3 



dQ i2dQ23dQ 31 S(^Q 3 - \lQl2+Ql3+Q 



{C3{Q3)-l) -~ (C2(Ql2)-l) - (C2(g23)-1) - (g2(Q3l)-l) 
^(C2(gi2)-l)(C2(Q23)-l)(C2(Q3l)-l) 



(6) 



r 



Note that the four terms in the numerator of the inte- 
grand in Eq. (|BJ involve four different (5-functions. This 
is clearly inconvenient when comparing to theoretical 
model calculations (for example those by Nakamura and 
Seki 0) which are not easily projected in the same way. 
More seriously, however, the mixing of different projec- 
tion algorithms for C3 and C2 in the numerator of Eq. © 
upsets the above-mentioned partial cancellation of the 
relative momentum dependences in the numerator and 
denominator of the integrand. It thereby acerbates the 
problem of introducing unwanted g-dependence into ra 
by projecting it. 

We demonstrate this by considering a very simple, 
spherically symmetric source characterized by the follow- 
ing (unnormalized) emission function: 



S{x,p) = Ep 6{t~tQ) e 



(7) 



Particles are emitted instantaneously at global time toj 
with a Gaussian momentum distribution which is inde- 
pendent of the emission point (no "x-p-correlations" ) . 
The factor Ep = ^Jrn^+p^ is introduced for calculational 
convenience so that it drops out from relations such 
as EpdN/cfip= J d'^xS{x,p) which relate the n-particle 
spectra Nn(p_u ■ ■ ■jJ'il) Eq. (0) to the emission hmc- 
tion S{x,p) [3, 0)13) 3- For simplicity we assume non- 
relativistic momenta, A < m. General arguments pre- 
sented in 0] suggest that for such a source the exact 



reduced three-particle correlator ra from Eq. Q should 
be equal to 2 and completely momentum-independent. 
The (nonrelativistic) two-particle exchange amplitude 
Pij — {d^ {pi)a{pj)) [2,Eitll!l for the source iQ reads 



S{x,K^j) 



(27ri?2)^e"i^ 



MEi-E,)to 



(8) 



where Kij = {pi+pj)/2 and Qij =Pi—Pj are the average 
and relative momentum of the pair ij. It allows to calcu- 
late the one-, two- and three-particle cross sections from 
the relations 0,0 

Ni{i) = pu, N2{ij) = piiPjj + Ipijp, 

^3(123) = 2Re(pi2P23/03l) + |Pl2pP33 + IP23IV1I 

+ \P3i\^ P22 + P11P22P33 (9) 
(where Ni{i) is shorthand for Ni{pi), etc.), yielding 



Ni{i)Ni{j) = (27ri?^)'^e~^e"3^, 



(10) 



and 



J 



Ari(l)7Vi(2)7Vi(3) 
iV3(123) 



(27ri?2)2 e~2^ 



(923 -931 ) 
e 12A2 



--R^(q23+q3i 



(11) 

(12) 



+e nxi— +e- "'"i2A2" +e- " 'sa^ 



912+923+912-923 



r 



After plugging these into equations (Q) and a bit of 
algebra confirms that the exact is indeed momentum- 
independent and equal to ra = 2. 

Let us now study the projected correlator r3(Q3) de- 
fined in Eq. ©. For non-relativistic momenta we can 
use —{Pi—PjY~{Pi~Pj)'^ = qij which allows to do the 
integrals in Eqs. (|4I5|) analytically. After transforming 
integration variables according to (Ppid^pj =d^Kijd^qij 
and d^pid^p2d^p3 — d^Kd^qijd^qjk (where {ij),{jk) in- 
dicate two of the three possible pairs (12), (23), (31), as 
suitable to exploit the cyclic symmetry of the three mid- 



dle terms in Eq. and further d^qijd^qjk — d^S^d^C, 

where C = Qij—Qjk and | = {Qij+Qjk)/'^, the integrations 
are straightforward, and we obtain for the projected cor- 
relators H4I5|I 



C2{Qij) 

C3{Q3) 



(13) 



2e 



— o S3 



- / dzz'Vl- 



2e 



-AS3 



6 e 



II (1^3) 



(14) 
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v^_p IE. erf {V^3 



(16) 



fa is plotted in Figure ^ as a function of S3 in the rele- 
vant range S3 ^ 3 where the two- and three-particle cor- 
relation functions C2{Qij) and C'3(Q3) show measurable 
deviations from unity [E S El- One observes a sig- 
nificant discrepancy to the exact result r3 = 2, and for 
large S3 one checks that fa grows exponentially nega- 
tive: f 3 ^ — I \Ji^eP^ /S3. At small relative momenta its 
leading (53-dependence is 



FIG. 1: (color online) The normalized true three-particle cor- 
relator ra (dashed, blue) and its projected version fz (solid, 
red) from Eq. ® for the simple source Q , as functions of the 
dimensionless variable E3 = Q\ (_R^-1/(4A^)) . 



Here z = cos0 arises from the angle between ^ and C,^ Ii 
is a modified Bessel function, and 



1 



4A2 



5^3 = ^12 + ^23 + S31 — Q\ \ — 



4A^ 



(15) 



The emission function S{x,p) is a Wigner density Q 
whose spatial and momentum widths must satisfy the 
uncertainty relation R^A^ > {h/2)'^. Hence, S.^ and S3 
are positive quantities. Inserting Eqs. I|13I14|I back into 
Eq. JBJ and introducing spherical polar coordinates, one 
angular integral can be done trivially while the radial in- 
tegral is executed with the help of the (5-function. The 
remaining integral over the polar angle yields an error 
function: 



^3(^3 



2 + 3 



£2 



i(S3/3) 



II (S3/3) 

^(^3/3) 



120' 



0(S3) 



(17) 



i.e. of 4 order in the Qij. The continued cancellation 
of 2"^^^ order terms is probably accidental and due to the 
particularly simple structure of the emission function Q. 
The coefficient in front of the S§ term is small, and the 
deviations of f3 from 2 become significant only at S3 >, 1. 
This may explain why Humanic did not see such a devia- 
tion in his Monte Carlo model study which explored 
fa in a restricted Q3 range and with limited statistical 
accuracy, even though he employed the same projection 
procedure 



IV. A DIFFERENT PROJECTION METHOD 

The (33-dependence in Eq. (|16|l is entirely artificial 
and generated by using different projection procedures 
in equations and (0) and mixing them in Eq. ((HJ. It 
can be easily avoided, however. Inserting the definitions 
^ into Eq. and multiplying both numerator and de- 
nominator by Ni{pi)Ni{p2)Ni{p3), rs can be brought 
into the form 



nno^ ^3(Pl,P2,P3) - N2{pi,P2)Ni{p3) - A^Z (P2 , P3)iVl (pi ) " (P3 , Pi )^1 (P2 ) + 27Vi (pi )7Vi (p2 ) A^l (P3 ) 

r3(123l = , 

V[^2(Pl,P2) - NMNi{p2)][N2{p2,P3) - m{p2)m{p3)][N2{p3,Pl) - iVl (P3)A^1 (Pl)] 
_ [NUM(pi,p2,P3)] 

[DEN(pi,p2,P3)] ■ ^ ' 

The numerator contains sums and differences of different kinds of triplet yields: the first term denotes real triplets 
where all three pions come from the same collision event, the last term mixed event triplets where all three pions come 
from different events, and the middle three terms subtract mixed-event triplets where two pions come from the same 
and the third from a different collision. For consistency with Eqs. (PJ, the distributions of both types of mixed-event 
triplets must be normalized such that they agree with the real triplets at large relative momenta where all quantum 
statistical correlations disappear: Eq. ||SJ) shows that for iy^j, pij vanishes in the limit \qij \ ^00. 

The denominator in (|18() involves a product of differences between yields of real pairs and mixed-event pairs (again 
normalized to the same total number of pairs), with pion momenta as indicated by the arguments. It can be rewritten 
as the yield of mixed-event triplets multiplied by a weight factor. 



[DEN(pi,p2,P3)] = A^l(Pl)A^l(P2)iVl(P3) ■ W(Pl,P2,P3), 



(19) 
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where the latter is computed from the (Coulomb corrected) two-particle correlation function C2{Pi,Pj) according to 



W(P1,P2,P3) = V[C2{Pl,P2) - 1][C2(P2,P3) " l][C2(p3,Pl) " 1] ■ (20) 

Fully six-dimensional analyses of the two-particle correlation function C2{Pi,Pj) = C2{Kij,qij) are available. It can 
be tabulated and interpolated, allowing for a straightforward calculation of the weight W(pi,P2,P3)- (Again the 
mixed-pair background must be normalized such that C2{\qij\~^oo) — 

We propose to construct directly from the ratio (|18|l instead of first constructing two- and three-particle correla- 
tion functions and then using |(2Jl. It is straightforward to project Eq. (|18l) onto the single relative momentum variable 
Qs (or onto any other desired combination of the three pion momenta) by binning the numerator and denominator 
separately in Q3 and taking the ratio: 

J d^Pld^P2d^P3 [NUM(pi,p2,P3)] S (Q3 - ^/-{Pl-P2? - {P2-P3Y - {ps^PlV) 

(^3)(Q3) - -4 ( ■■ (21) 

/ d3pid^p2d^p3 [DEN(pi,p2,P3)] S [Qs - V-(.Pl-P2)^ - {P2-P3r - {PZ-pif) 



Note that, by symmetry, the three middle terms in the 
numerator of Eq. H18|) give identical contributions to 
Eq. H21|) . The numerator is thus obtained by first run- 
ning separately through all real triplets, all triplets with 
two particles from the same event and the third parti- 
cle from a different event, and all mixed-event triplets, 
binning each of these in Q^. After normalizing the bin 
contents of the mixed-event triplets by the appropriate 
factor which ensures that their distributions agree with 
that of the real triplets at large \qij\ (note that this does 
not imply that they necessarily agree at large (53!), one 
then subtracts, for each value of Qa, from the bin con- 
taining the real triplets three times the bin content from 
the half-mixed triplets and adds twice the content from 
the fully mixed triplets. For the denominator one runs 
over all fully mixed triplets, calculates for each of them 
the weight W and bins these weights in Q^. At the end 
one normalizes the bin contents by the same factor as for 
the fully mixed triplets in the numerator. 

While this projection procedure does not, of course, 
completely sidestep the unavoidable interference with 
the cancellation of the leading relative momentum de- 
pendences between numerator and denominator, at least 
it preserves the momentum-independence of in the 
"trivial" case iQ of a spherically symmetric, static 
source without x-p-correlations: For the source |(7J), 

(r-3)(Pl,P2,P3) = ?'3(P1,P2,P3) = 2. 

A technical complication for the data analysis arises 
from the repulsive final state Coulomb interaction be- 
tween the emitted charged pions which must be corrected 
for before computing the reduced three-pion correlator 
r^. This problem poses itself similarly for both projec- 
tion techniques and can be dealt with in the same way in 
either approach, by using Coulomb correction factors ex- 
tracted from a self-consistent Coulomb correction of the 
two-pion correlation functions. We refer to Refs. 0, Q 
for details. Note that in the form given in Eq. (|18|l . the 
Coulomb correction must be applied to the real pairs 
and triplets, and not to the mixed-event ones. 



V. CONCLUSIONS 

We have shown that, for a simple spherically sym- 
metric, non-expanding source with sudden particle emis- 
sion, the standard projection procedure, as up to now 
applied in the experimental analysis of three-pion corre- 
lations, introduces an artificial momentum dependence 
into the normalized true three-pion correlator r^, even 
though the exact result for is completely momentum- 
independent. This suggests that the observed depen- 
dence of fa seen by WA98 in Pb-|-Pb collisions at the SPS 
111 and by STAR in Au-hAu collisions at RHIC 6] may 
be significantly affected by artificial projection effects. 
Even though projection-induced relative momentum de- 
pendence of ra can in general not be completely avoided, 
it can be minimized by using the new projection proce- 
dure suggested in this paper. This procedure computes 
ra in a more direct way and applies the same binning to 
all terms. It preserves the momentum-independence of 
for a symmetric non-expanding source. The new projec- 
tion algorithm can also be easily employed in theoretical 
model calculations, rendering the comparison with data 
(but not necessarily the interpretation of the extracted 
source parameters) more straightforward. 
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